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Abstract 


The paper is devoted to the physical parameterization of the top-of-atmosphere reflection function. The ac- 
curacy of the parameterization is checked against exact radiative transfer calculations in a cloudy atmosphere 
for various cloud-top-heights, cloud optical and geometrical thicknesses, solar illumination and surface reflec- 
tion conditions. It was found that the error of approximation is smaller than 5% for most cases studied at the 
wavelength interval 758-768 nm, which corresponds to the oxygen A-band. This band is routinely used in 
cloud-top-height retrievals. The model proposed can be used for cloud-top-height and cloud geometrical thick- 
ness retrievals. This allows to avoid a standard look-up-table retrieval scheme, involving complex numerical 
procedures. 
© 2003 Elsevier Ltd. All rights reserved. 
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1. Introduction 


The radiative transfer in the Earth-atmosphere system is usually studied numerically, solving the 
well-known integro-differential equation [1]. This is mostly due to the necessity to account for 
multiple light scattering by aerosol particles and molecules. 

The appearance of clouds generally complicates the situation due to the necessity to consider the 
random structure of cloudiness [2]. On the other hand, the assumption of an optically thick homoge- 
neous plane-parallel cloud layer, which is of importance in most of practical situations [3], leads to 
the simplification of the problem. This is mostly due to the fact that highly variable reflectance of 
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light from the surface and boundary-layer aerosol beneath the cloud plays only a minor role in the 
total solar light reflectance from cloudy scenes. Thus, they can be accounted for using rather coarse 
simplifying assumptions and models. 

The task of this paper is to propose a simple semi-analytical model for the calculation of the 
top-of-atmosphere (TOA) reflectance of a underlying surface-atmosphere system, accounting both 
for aerosol and cloud scattering. The atmospheric gas absorption is also taken into account. Main 
attention is given to the calculations for the spectral interval, determined by the oxygen A-band 
(758-768 nm). This band is important for the cloud-top-height determination from spaceborne in- 
struments. Indeed, the reflected light intensity is extremely small in the center of this band for 
cloudless conditions. This is due to the fact that photons can penetrate down to the underlying sur- 
face and be absorbed by oxygen molecules on their way to the surface and back to a receiver. For 
cloudy atmospheres the situation is different, however. Imagine, that cloud is present at the height 
10 km above the ground level. Then most of photons are reflected by the cloud and do not penetrate 
in the layer underneath the cloud to be absorbed there. This will lead to an increase of reflected light 
intensity in the center of oxygen A-band as compared to the case of cloudless atmosphere. Clearly, 
the depth of this line J depends on the cloud-top-height. The value of J is larger for clouds having 
smaller cloud-top-heights. 

Note that other absorption bands can be employed for the same purpose (e.g., oxygen B-band, O4 
absorption bands, etc.). Our model can be also used to treat all these different cases (both in and 
out of gaseous absorption bands). 

Moreover, the simplified solutions obtained are of importance to clarify the physical processes, 
which contribute to the TOA reflectance. They can be used also for the development of analytical 
and semi-analytical cloud retrieval schemes. 


2. Theory 


The schematic representation of the geometry considered is given in Fig. 1. Unpolarized solar light 
illuminates the aerosol-gaseous atmosphere in the direction Qo(Wo, po) and a receiver detects the 


reflected light intensity Z in the direction, specified by the angle Q(v,@). Here Vo and V are zenith 
incident and observation angles respectively, pọ and @ (not shown in Fig. 1) are correspondent 
azimuth angles. We will assume that @) = 0. Then the azimuth difference between incident and 
reflected beams directions is given by the angle ọ. 

A cloud layer with top height H and geometrical thickness AH (see Fig. 1) is introduced in the 
aerosol-gaseous atmosphere. Both AH and H can be in principle arbitrary. However, we assume 
that AH > 51, where / is the photon free path length in a cloud. For a typical value / = 20 m, it 
follows: AH > 100 m. 

It means that only cases of optically thick clouds with the optical thickness te = AH/] >5 
are considered. Note, that thinner clouds are highly inhomogeneous both in vertical and horizontal 
directions. Therefore, the model of a plane-parallel homogeneous cloud layer is characterized by a 
large uncertainty for thinner clouds. We avoid these difficulties, assuming from the very beginning 
that te > 5. The surface is assumed to be Lambertian with albedo A. 

The TOA reflectance Ra,(vo, V, p), measured by a radiometer or spectrometer with a spectral 
resolution A4, is related to the reflected monochromatic light intensity (A, Vo, V, p) by the following 
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Fig. 1. The geometry of the problem. 


relationship [1]: 


Isa, V0, 
Rað, Vo, 0) = a 


Fa, cos vo ° 


(1) 


where 


2+42/2 2+44/2 
Itobo SARA oat, fOO SONA, 
2—41/2 A—Ai/2 

where f(A, 2’) is the radiometer spectral response function and zF',, cos Yo is the incident solar flux 
density on the TOA in the spectral range AZ. Value (1) is routinely measured by various radiometers 
and spectrometers, placed on satellite platforms [4]. 

The shape and half-width of the function f(A, 1’) is only of a minor importance for measurements 
outside of atmospheric gases absorption bands. This is not the case, of course, for absorption bands of 
atmospheric gases (e.g., O2, H20, CO2, N20, CH4), which have an extremely complicated structure. 
The task of this paper is to model TOA spectra both outside and inside gaseous spectral bands. Thus, 
the precise shape of f(A, A’) cannot be ignored. We, therefore, develop here the approximation for the 
monochromatic intensity /(/, Vo, V, p) or, more precisely, for the monochromatic reflection function 
R(A, 0,90, 9) = 104, 0, V0, 9)/F cos Vo. This allows to use the results presented below for any types 
of the response functions f(A, 2’). 

The TOA reflectance R can be presented as a sum of two terms: 


R= Ri + Tk), (2) 


where R describes light scattering and radiative transfer in the atmosphere above a cloud (with 
account for both gaseous and particle scattering and absorption). The value of R is due to 
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cloud-underlying atmosphere and surface contribution (see Fig. 1). We omit arguments in func- 
tions in Eq. (2) for the sake of simplicity. The multiplier T accounts for the extinction of the direct 
solar light on the path from the sun to the cloud top and from there to a receiver. 

The scattering of light above clouds is rather weak. Thus, we will calculate the value of R, in 
the single scattering approximation [1]: 


1 h—H 7 7 
Rae | dz(o%alZ)pP® (0) + assalz) p^(0)) exp — TEET! + n7’), (3) 
0 
where o®:A(z) are Rayleigh and aerosol scattering coefficients, respectively, p®4(0) are Rayleigh and 


aerosol phase functions, t is the optical thickness of atmosphere above cloud along the vertical axis 
OZ (see Fig. 1), directed to the Earth surface, and h is the TOA height, assumed to be equal 60 km 
in this study. The value of t includes the contribution from the molecular and aerosol scattering and 
absorption. Also we introduced here the values of € and y : č = cos Vo, n = |cos ù|. 

It should be pointed out that coupling between atmospheric layers above and below cloud-top-height 
is neglected in Eq. (2). To account partially for this coupling and multiple light scattering above 
cloud top, we assume: 


T = exp[ — T'(&7! + 47')], (4) 
where 
N h-H 7 å 
= Cy. (z)e; . 5 
T >, I OTOLE (5) 


Here c$ C) is the ith gas absorption cross section, N the total number of gases present and c; the 
ith gas concentration. Only extinction due to light absorption by gases is present in 1’. 

Actually, the value of t’ should include also extinction by aerosol particles (the optical thickness 
tê) and molecules (the optical thickness tÈ). We neglect them in Eq. (5). This results in increase of 
T and the term TR, in Eq. (2), which partially compensates for multiple light scattering in the layer 
above cloud, which otherwise is completely neglected in this study. Note, that this is a standard 
method for a partial account of multiple light scattering effects in the problems of the atmospheric 
correction. However, the definition of t’ in Eq. (4) differs, depending on the problem and spectral 
range studied [5-7]. 

Now we consider the term Rz. The general expression for this term for the cloud-underlying 
Lambertian surface is well-known [4]. It has the following form: 


At, (Vo )te(P) 


Rv, Vo, o) = Rv, Vos o) + l — Ar, > 


(6) 
where R.(o, V, p) is the cloud reflection at Æ = 0, ¢,(v) is the diffused transmittance of a cloud 
layer [4] and rẹ is the spherical albedo of a cloud. We also define the cloud total transmittance: 
t-=1—Pre. 

Convenient approximate equations for Re, te, and re were obtained in many studies [1,4,8—11]. 
Those given in [11] and used in this study are summarized in the appendix. Eqs. (2), (3), (6) do 
not completely solve the problem we face. We have an aerosol—gaseous medium between underlying 
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surface and cloud (see Fig. 1). It will influence relationship (6), obtained for the case of a transparent 
layer between cloud and a ground surface. 
To approximately account for this influence we substitute 4 in Eq. (6) by the following equation: 
At? 
1 — Ar,’ 


A* =r + (7) 
where rą is the spherical albedo of aerosol layer between cloud and a underlying surface and ta 
is the aerosol total transmittance (ta = 1 — ra for non-absorbing aerosols). Eq. (7) is analogous to 
Eq. (6). To show this, we integrate Eq. (6) over angles Vo, V and g. Then it follows [10]: 


At? 
1 — Arg 


n =r: + (8) 
which is similar to empirical Eq. (7). Note that r2 =r. (4=0). In essence, approximation (7) allows 
to substitute an aerosol-underlying surface layer by an effective Lambertian surface with albedo A*. 

Clearly, we have at ra = 0 : A* = A, where it is assumed that most of atmospheric aerosols are 
weakly absorbing and tą = 1 — ra. Thus, for calculation of A* we should find the value of ra (see 
Eq. (7)). Accounting for the fact that aerosol optical thickness is usually small and contribution to 
the TOA reflectance from the cloud underneath is low, we expect that rather coarse approximations 
for the calculation of 7, can be used. In particular, one can obtain in the single scattering approxi- 
mation [12]: 


1 1 
— d dy F(é,n) p(é,n), 9 
, on | eae f nan FEM mE (9) 
where 
1 expl—t8(€! 49] 
F(En)= = (10) 
and 
1 2n 
PEM f Odo. (11) 
Here 
0 =arccos(—én + \/(1 — €2)(1 — 7?) cos o) (12) 


is the scattering angle and tê is the aerosol optical thickness beneath the cloud. 
It was shown [12] that under assumption tê — 0 Eq. (9) is reduced to a simple function expressed 
via the elliptic integral 


i a 13 
K(p)= ——. 
m | V1 -— psin’ 0 ma 


Namely, it follows: 


Fa = 20 fr", (14) 
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where 

z_ l-g { 2K(g°) 

= 1 15 

b=- { ~ (15) 

is so-called backscattered fraction and g is the asymmetry parameter: 
1 T 
v=; / p(@) sin 0 cos 0 dé. (16) 
0 


Eq. (15) is exactly applicable to the case of the Heney—Greenstein phase function: 


1-9? 
(1 + g? — 2g cos 6)3/2 


p(0) = = X_ (21 + 1)g'Pi(cos 0), (17) 
l=1 


where P;(cos 0) is the Legendre polynomial. However, it can be used with a high accuracy for other 
phase functions if g is known [12]. 

We present integral ra (9) in Fig. 2 as the function of t^ for different values of g. The value of 
ra (tå =2) approximately coincides with ra (t^ — 00) for the range of g studied. So ra (t^ = 2) 
gives the contribution of the single scattering to a light reflection from a semi-infinite medium. We 
also see that the linear approximation (14) is valid for values of tê =0.05—0.1, depending on g. For 
larger values of tê the dependence r(t) becomes non-linear. In particular, we found approximating 
the function r,(t*), given by Eq. (9), at g = 0.77 and tê < 0.3 (see Fig. 3): 


ra(T) = @o[0.265t* — 0.865(t*)? + 1.164(c4)°]. (18) 


To complete the model, we should account for the light absorption of diffused light fluxes by 
gases present in a layer between a cloud and an “effective” underlying surface (see Eq. (7)). To do 


spherical albedo 


06 08 10 12 14 16 #18 20 
optical thickness 


Fig. 2. The dependence r,(t*), obtained with Eq. (8) at g = 0.55, 0.66, 0.77 and 0.88. 
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Fig. 3. The dependence r,(t*), obtained with Eq. (18) (line) and Eq. (8) (starts) at g = 0.77. 
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Fig. 4. The relative error of Eq. (20) as compared to Eq. (19) as function of tc”. 


so we multiply te(Vo) and ¢.(v7) in Eq. (6) by the flux transmittance [13,14] 
1 * 
y=2 | aččexp -5 (19) 
0 č 


where t* is the optical thickness due to absorbing gases in a layer between a cloud bottom and a 
ground surface. Eq. (19) can be substituted by 


y = exp(— Tr"), (20) 


where Y ~ 1.7 is the so-called diffusivity factor commonly used in the transmission modeling. The 
accuracy of Eq. (20) is better then 2% at t* < 0.4 (see Fig. 4). For larger values of t* one should 
perform a simple integration (19) (if the high accuracy is needed). Note that y — 0 as t* — oo. 
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Summing up, we have the following final approximate relationship for the TOA reflectance: 


A*t (Ot (V 
R(9, vo, P) = Ri (9, 80, P) + |Re(9, Vo, p) 4 a Í 7, (21) 


where Ri, 7,A*,y are given by Eqs. (3), (5), (7), (20) and approximate analytical formulae for 
Re, fc, tc(Vo) are given in the appendix. 

The importance of Eq. (21) is not only in the possibility to substitute the full radiative transfer 
calculations in the Earth—cloudy atmosphere system by an approximate formula. Most important 
is that it offers a simple way for the semi-analytical solution of the inverse problem (e.g., the 
cloud-top-height determination). This will be a subject to our next publication, however. 


3. The forward model 
3.1. The radiative transfer model 


Owing to the phenomenological derivation of Eq. (21), it is not possible to derive the errors intro- 
duced, using only the model presented above. Clearly, for this we need to run a full radiative transfer 
code for the Earth atmosphere-underlying surface system with realistic distributions of atmospheric 
gases and aerosols. For this study we have used the SCIATRAN radiative transfer code, developed 
at Bremen University [15]. The SCIATRAN code originates from the well-known and highly used in 
trace gas retrievals GOMETRAN model [16]. The test of SCIATRAN accuracy can be found in [17]. 
In particular, it was found using the DISORT model [18], that the accuracy of SCIATRAN is better 
than 1% for most of geometries and any cloud optical thickness. Note that calculations according to 
the Mie theory [10] are used as an input to SCIATRAN radiative transfer code. Those allow to find 
phase function, single scattering albedo and extinction coefficient as functions of particles’ refractive 
index and particle size distributions. The cases of non-spherical particles can be also treated, but 
then the Mie theory should be substituted by light scattering theories valid for such scatterers [10]. 

We found that the standard version of SCIATRAN [16,17] has a poor accuracy (errors up to 
5%) for water clouds at the rainbow and glory scattering regions. This is due to a limit number of 
Legendre polynomials, used in the cloud phase function representation. To overcome this problem, 
the following correction of SCIATRAN [16,17] was developed. First of all we have calculated the 
contribution of the singly scattered light /,,, using the SCIATRAN model. Then we have subtracted 
Iss from the result obtained and added the exact value of /;,, which has a simple analytical form 
[10]. We found by comparison with published results [19] that the accuracy of this new version of 
SCIATRAN for cloudy atmospheres is better than 1% also for rainbow and glory geometries. 

We restricted the comparison by an extremely narrow spectral range of oxygen A-absorption band 
(758-768 nm). It was done by purpose. This spectral range is usually applied to the cloud-top-height 
determination [20-23]. Thus, Eq. (21), proven to be valid, can be applied for the same purpose. 
Also the value of the reflection function and the atmospheric transmittance is changed practically 
from 0 to 1 in the band and its vicinity. This allows us to consider a very broad range of TOA 
reflectances, even if consider only this small spectral interval (AA = 10 nm). 

The final note is related to the comparisons of exact results with approximation (21). Namely, we 
do not present here the comparisons for the case of monochromatic radiances. This is mainly of an 
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Fig. 5. The dependence o%,(z) (broken line), o%,(z) (solid line), used in calculations (the wavelength 2 = 760 nm). 


academic interest. Instead, we have used Eq. (1) to find Ra, at AA=0.35 nm. This corresponds to the 
spectrometer GOME [26]. Note, that other spectrometers in current operation have larger values of 
A2. We found that the accuracy of Eq. (21) increases with AJ. Thus, all estimations presented here 
can be considered as upper error estimations (at least as currently flying spectrometers are concerned, 
including most advanced SCIAMACHY spectrometer on board of ENVISAT, having A2 = 0.48 nm 
at 2 = 604-805 nm [24]). The response function has been taken in the following form: 
ql 4 4 4 
Gane 1/Aa, X e[A— Ad/2,A+ Ad/2], 


0, otherwise, 


which is a simple box function. 
3.2. The model of atmosphere 


The vertical distributions of aerosol and Rayleigh scattering coefficients through terrestrial atmo- 
sphere used in this study is shown in Fig. 5 at the wavelength 760 nm. The values of tô and tÈ were, 
respectively, 0.26 and 0.027 at this wavelength. Aerosol properties are taken from the LOWTRAN 
model [25]. The tropospheric aerosol model used corresponds to the case of maritime winter aerosol 
with ground visibility 23 km and boundary layer humidity 80% [25]. The stratospheric aerosol model 
corresponds to the so-called background aerosol [25]. 

We clearly see the enhanced aerosol layer around 20 km and highly pronounced boundary aerosol 
layer in the range 0-5 km above from the ground in Fig. 5. The vertical structure of the gas absorp- 
tion coefficient is shown in Fig. 6 for wavelengths 758.025, 761.525 and 763.525 nm. The correspon- 
dent O2 optical thicknesses are 0.0029, 43.1 and 0.26 for these wavelengths. Note a considerable 
change of the optical thickness over a small spectral interval (5.5 nm). 

The strong absorption (see broken lines in Fig. 6) is due to oxygen. Solid line presents compar- 
atively weak absorption by other gases (mostly ozone with the optical thickness equal to 0.004), 
which is identified by the peak around 20 km (see Fig. 6). The data presented in Fig. 6 and also at 
other wavelengths used in calculations (not shown explicitly in Fig. 6) are calculated utilizing the 
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Fig. 6. The dependence of gaseous absorption coefficient on the altitude for selected values of 2 = 758.025, 761.525 and 
763.525 nm. 


HITRAN 96 database of line parameters [26]. We assumed that the lines shapes are given by the 
Voigt function. All calculations are performed for the case of temperature and pressure profiles for 
the 55° latitude (January) according to the model [27]. 

The phase function for aerosol was taken in form (17) with values of g, depending on the 
height in a discrete way. In particular, we assumed that g = 0.7731 at z* < 2 km, g = 0.6614 at 
z* € [2 km, 10 km], g=0.6448 at z* € [10,30 km], g = 0.6868 at z* > 30 km. Here z* is the attitude 
above the sea level. 

The Rayleigh phase function was assumed to be constant with height and given by a familiar 
expression [28]: 


r= 


3 2 
rina a p) cos’ 0], (22) 
where p is the depolarization factor. 

The aerosol is assumed to be weakly absorbing with absorption coefficient given by the Lowtran 
model [25]. We do not present the correspondent dependence here due to its small influence on the 
TOA reflectance in the framework of the model under study. 

The results presented below were obtained for the Deirmendjian’s Cloud C.1 droplet distribution 
f(a) = Aa®e~ ©, where A = const and aj) = 4 um. We have used 50 moments of the expansion 
of the phase function on the Legendre polynomials. This (combined with the ô-M approach [1]) 
allows accurately account for the phase function angular dependence in our radiative transfer code. 

The vertical distribution of the atmospheric single scattering albedo used in calculations is given 
in Fig. 7 for wavelengths 758.025, 761.525 and 763.525 nm. The broken line with dots near the right 
boundary of Fig. 7 indicates the single scattering albedo, which is due to aerosol absorption only (in 
the assumption that os. = 0). We see that the aerosol absorption is indeed small and w  ~ 1.0 for 
the aerosol medium. The solid line presents the single scattering albedo outside the oxygen A-band. 
It is mostly due to ozone absorption. We see that wo > 0.9 at this wavelength. It is close to one 
in the region of heights 4-5 km, where cloud is placed. This is due to non-additivity property of 
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Fig. 7. The dependence wo(z) for three wavelengths, used in this study. The curve, with highest values of wo(z) corresponds 
to the artificial cal o$, = 0. 


the single scattering albedo, combined with the large values of scattering coefficients of a cloudy 
medium. 

In comparisons presented below the values of H and AH (see Fig. 1), however, are not fixed, but 
varied. It allows us to study the influence of those variations on the accuracy of main 
Eq. (21). Also we have varied the values of A,Vo, cloud optical thickness, the aerosol type and 
cloud microphysical parameters. The zenith observation angle Ý was fixed and assumed to be equal 
0° (nadir observation). 

Summing up, we underline that our atmospheric model is very close to a realistic scenario in the 
Earth atmosphere. And yet we found that the accuracy of Eq. (21) is better than 5% for most of 
cases studied. This is close to typical measurement errors. Some results of comparisons performed 
are given below. 


4. The accuracy of the approximate model 
4.1. The cloud optical thickness variation 


First of all we study the accuracy of our approximate model in respect to the cloud optical 
thickness variation (see Fig. 8). The results obtained are also of a general interest. 

The calculations in Fig. 8 were performed for values of cloud optical thickness 5—100, which 
covers the most frequent range of cloud optical thickness change [3]. We see that the TOA reflection 
function is extremely sensitive to the value of cloud optical thickness outside the oxygen absorption 
band, located at 4 = 758.5 nm. It varies in the range 0.25—0.8 for the case studied. This feature is 
used in the cloud optical thickness determination [29,30]. 

The dependence of the TOA reflection function on the cloud optical thickness inside the oxygen 
absorption band is also of importance. In particular, for the case studied here the TOA reflection 
function in the center of O) A-band increases from 0.05 to 0.1, while cloud optical thickness changes 
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Fig. 8. The dependence of the TOA reflection function R(0°,60°) on the wavelength in the oxygen A-band for the 
atmospheric model, specified in the text at AH = 1 km, H = 5.5 km, A =0 and several values of cloud optical thickness 
Te = 5, 10,21, 100 (1—exact calculations for the Earth-cloudy atmosphere system, 2—calculations according to Eq. (21)). 
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Fig. 9. The relative error of Eq. (21) ¢ = 100(1 — R/R*), where R is obtained from Eq. (21) and R* from exact radiative 
transfer calculations for the case given in Fig. 8. 


from 5 to 20. Thus, to derive cloud-top-height from the depth of the absorption line around 761 nm, 
one should know also the value of te. 

We see that approximate Eq. (21) describes correctly all characteristic features of the spectral 
behavior of the TOA reflectance function. The error of Eq. (21) is given in Fig. 9. 

It follows from Fig. 9 that the accuracy of Eq. (21) is in the range [ — 5%,3%] for the cloud 
layer of 1 km thickness and the top height equal to 5.5 km. Such a high accuracy is hardly expected, 
taking into account the phenomenological derivation given above. It, therefore, confirms our physical 
arguments concerning the significance of different contributions to the TOA reflectance. 
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Fig. 10. The same case as in Fig. 8, but at 4=0, AH =0.25 km, te=5 and several values of H=0.5, 3.5, 5.5, 7.5, 9.5, 11.5 km. 


Similar accuracy was obtained for other solar angles in the range 0-70°, other cloud droplets 
average radii in the range 4-16 um and other realistic extinction coefficients of a boundary aerosol 
layer with the meteorological visibility range 10—50 km (not shown here). 


4.2. The cloud-top-height variation 


The dependence of R on the cloud-top-height is given in Fig. 10. The results were obtained at 
t=5, AH =0.25 km, vo = 60°. We see that the TOA reflectance increases with cloud-top-height from 
R=0.01 (at H = 0.5 km) till R=0.14 (at H = 11.5 km) in the center of the absorption line. The 
increase is in more than order of magnitude. This fact is used for cloud-top-height monitoring from 
satellites [20—23]. 

The physical origin of the TOA reflectance dependence on the cloud type height around 4=761 nm 
is clear. It is due to high reflection of light from a cloudy medium. This does not permit photons 
to penetrate to deep atmospheric layers and be trapped there by oxygen molecules. 

The results of comparisons of calculation of the TOA reflection function, using Eq. (21) and 
exact theory, is presented in Fig. 10. We see that Eq. (21) correctly describes the dependence of R 
on the cloud-top-height. The correspondent errors are given in Fig. 11. Note that errors are smaller 
than 4% even for low clouds with cloud-top-height equal to 0.5 km. The error decreases down to 
2.5% for clouds with cloud-top-height 11.5 km. It suggests that the single scattering approximation 
is a valid approximation even if clouds are absent, A = 0, and À = 758-768 nm. Most probably this 
is due to a high light absorption in oxygen A-band, which reduces multiple scattering effects in a 
clear atmosphere. This is confirmed by numerical experiments at 1 = 758-768 nm and A = 0. Note, 
that the accuracy of single scattering approximation decreases with the value of A due to additional 
illumination of atmosphere from the ground. 

In conclusion, we underline that the accuracy of Eq. (21) inside the oxygen band depends also 
on the cloud geometrical thickness. It follows from Fig. 11 that the accuracy of Eq. (21) is 0.5% 
in the center of the oxygen band at H = 5.5 km. This is in contrast with 5% error for the same 
conditions (e.g., the same value of te), but different cloud geometrical thickness, given in Fig. 9. 
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Fig. 11. The relative error ¢ for the case presented in Fig. 10. 


Cloud geometrical thickness for data given in Fig. 9 is 1 km as compared to only 0.25 km for 
the case, given in Fig. 11. It means that the cloud extinction coefficient is large for the case, 
given in Fig. 11. Absorption coefficient is, however, the same for both cases. It means that the ratio 
absorption/extinction is larger for the case given in Fig. 9. This is a genuine reason for the difference 
in errors for both cases. Eq. (21) becomes more accurate with decrease in the level of absorption. 
This is mostly due to the application of formulae given in the appendix. Their accuracy decreases 
with 1 — wo, where Wo is the single scattering albedo [11]. 

Data given in Fig. 10 can be used for the cloud-top-height estimation from measurements of the 
reflection function. Note that it is of advantage to use relative values in the estimation of H. In 
particular, we can introduce the relative band depth: 


R(A1) — R(%2) 


= 2 
- R(A\) (23) 
or 
@=1-Q", (24) 
where 
R(A1) 
2 
RO2) (25) 


Here 4; is the wavelength outside of the band and 22 is the wavelength inside the absorption 
band. Clearly, we have: Q > 1. The dependence of Q on the cloud-top-height is given in Fig. 12a 
both using SCIATRAN radiative transfer solver and approximate equation (21). We assumed that 
2; = 758 nm and A, = 763 nm. Also we used: AH = 0.25 km(t = 5) and AH = 1.0 km(t = 20). 
Other input parameters coincide with those used in Fig. 8. We found that the relative error of Eq. 
(21) is less than 3.5% for the case given in Fig. 12a both for t equal to 5 and 20. This is very 
promising result as far as cloud-top-height retrievals, using our approximation, are concerned. They 
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Fig. 12. (a) The dependence of the ratio R(758 nm)/R(763 nm) on cloud-top-height H at t = 5(AH = 0.25 km) and 
t=20(AH =1 km). Other conditions as in Fig. 8. Lines give approximate results. Points are obtained using SCIATRAN. 
(b) The dependence of the absolute cloud-top-height error Y on the value of the cloud-top-height H obtained with 
Eq. (26) (lines 2,4) and Eq. (28) (lines 1,3). Other input parameters as in (a). 


will be treated in a separate publication. However, for the sake of completeness, we show in Fig. 
12b the absolute accuracy of the cloud-top-height determination with approximate theory developed 
here, using SCIATRAN data presented in Fig. 12a as measured ratios. Solid lines in Fig. 12b give 
absolute cloud-top-height error ¥ = Hı — Hy, where H, is the measured (or, more exactly, obtained 
from the exact SCIATRAN model) cloud-top-height and H, is the cloud-top-height obtained from 
the application of our approximate theory. More precisely, in calculations of e, we have used the 
following relationship: 


(26) 


where @ is derivative of @ with respect to H, 0, is the value of © obtained from the SCIATRAN 
(see Fig. 12a) and @, is the same value, which follows from our approximate theory (see Fig. 12a). 
We see that the absolute error is negative. Therefore, our approximate theory overestimates the value 
of H. This overestimation is not larger than 0.25 km for cloud-top-heights in the range 0.5—11.5 km. 
Clearly, the correspondent relative error decreases rapidly with H. 

Note that the value of ¥ can be reduced if more spectral channels (and not only 2,22) are used. 
To illustrate this we introduce the logarithm 


a} 
Rx(A) J’ 


E(A) =1n { (27) 


where subscripts have the same meaning as in Eq. (26). Then the value of Y can be obtained 


minimizing the quadratic form 


din Ro(A) y 
dH 


2 


F(A) > min. (28) 
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Fig. 13. The same as in Fig. 8, but at te = 20 and AH = 1 and 5 km. 


The results of this procedure at the wavelength range 758-768 nm (other parameters as in Fig. 12a) 
are given in Fig. 12b by dashed lines. We see that this approach allows to reduce the absolute error 
of the cloud-top-height determination to values smaller than about 150 m for values of H in the 
range 0.5-12 km. 


4.3. The cloud geometrical thickness variation 


The cloud optical thickness Te is defined as the product of the cloud geometrical thickness AH 
and the cloud extinction coefficient Gext. It means that clouds, having the same value of te, but 
different values of AH, will have different values of the extinction coefficient, as discussed above. 
In particular, the extinction coefficient is larger for geometrically thinner clouds at te = const. As 
it was discussed above, it leads to smaller values of the single scattering albedo for geometrically 
thicker clouds. Thus, we note that the increase in cloud geometrical thickness at te = const should 
lead to the decrease of the TOA reflection function. 

These arguments are confirmed by data given in Fig. 13, obtained for various values of AH 
and t = 20,09 = 60°, A = 0. We see, therefore, that the dependence R(AH) cannot be ignored. For 
instance, the increase of AH from 1 to 5 km leads to a two-fold decrease of the TOA reflection 
function (from 0.1 to 0.05 at 4 = 760.5 nm). Large variations of R occur at the wavelength 763 nm, 
where they reach 30%. Thus, the cloud-top-height retrieval algorithms should incorporate also the 
value of AH estimation (e.g., from the water vapor absorption band around 938 nm [23]). On the 
other hand, knowing the value of H (e.g., from the Ring effect [31]), one can obtain also AH from 
measurements in the oxygen A-band. The accuracy of Eq. (21), given in Fig. 14, appears to be 
high as well for this case. It is in the range [ — 3%,3%] for AH = 1 and 3 km. The error sharply 
increases around 4 = 761 nm and AH = 5 km. However, this case is presented for illustration pur- 
poses only. First of all clouds with t = 20 do not have the thickness around 5 km. They are much 
thinner. 
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Fig. 14. The relative error ¢ for the case presented in Fig. 13. 
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Fig. 15. The same as in Fig. 8, but at AH = 0.25 km, H = 9.5 km, te = 5 and several values of A = 0, 0.2, 0.4, 0.6, 0.8. 
4.4. The surface albedo variation 


We present results of the calculation of the reflection function in the oxygen A-band for different 
Lambertian ground albedo and t = 5 in Fig. 15. We assumed that cloud geometrical thickness is 
equal to 0.25 km the cloud-top-height is equal to 9.5 km and Vo = 60°. Note, that effects of ground 
albedo are of smaller importance for larger values of cloud optical thickness. 

In the center of the oxygen A-band, the variation of the surface albedo almost does not influence 
the TOA reflectance. Outside the oxygen absorption band the influence of the surface albedo is 
important. We found that the variation of albedo from 0 to 0.8 leads to approximately the same 
variation of TOA reflectance as due to the variation of the cloud optical thickness from 5 to 100 
(compare Figs. 8 and 15). 

The comparison of exact and approximate calculations show that Eq. (21) grasps all important 
features of TOA reflectance variation inside the oxygen A-band, depending on the ground albedo 
(see Fig. 15). 


52 A.A. Kokhanovsky, V.V. Rozanov! Journal of Quantitative Spectroscopy & Radiative Transfer 85 (2004) 35-55 


error, percent 


758 760 762 764 766 768 
wavelength, nm 


Fig. 16. The relative error ¢ for the case presented in Fig. 15. 


The accuracy of Eq. (21) is smaller than 8% (see Fig. 16) even for very bright surfaces with 
A=0.8, which roughly corresponds to snow fields. It is in the range 2% for A =0.4, which is typical 
for deserts. 


5. Conclusion 


In conclusion, we present here a simple physical parametrization (21) for the TOA reflectance. 
This parametrization can be used to avoid calculations of the TOA reflectance over cloud scenes, 
using exact radiative transfer codes. 

The error of Eq. (21) in most cases is smaller than 5%, which allows its use also for the solution 
of the inverse problem, namely for the determination of t, H and AH. 

Most of comparisons have been given for a practically important spectral band 758-765 nm. 
However, there is no reason to expect that model will fail to reproduce experimental results for 
other wavelengths in visible and near-infrared (at least up to 2 = 2.4 um). 

Both approximate solution and exact calculations were used to study the dependence of the TOA 
reflection function on various parameters, such as cloud-top-height, cloud geometrical and optical 
thickness, ground surface albedo, illumination conditions, etc. In particular, we found that cloud ge- 
ometrical thickness value cannot be ignored in cloud-top-height retrievals, using the oxygen A-band. 
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Appendix A. The cloud reflection function 


Reflection function of a homogeneous optically thick (optical thickness te > 5) cloud over a black 
surface (e.g., ocean) is given by the following analytical solution [8]: 


RO, vo, P) = RY(V, Vo, P) exp(—uy(1 — 0.05 y)) — (t exp(—x) — t“) exp(—y)K (Yo) K (0), 
(A.1) 
where 
inh 
t= eee ee t* = (4.86 — 13.08 cos Jy cos 3 + 12.76 cos” 9o cos” 9)/t3, 
sinh[ y(x + z)] 


K(0) = 3(1 + 2cos¥), 
1- Wo 


a= 1.07, x=yz, y=4 30 = 9)’ 2= 51 gro 


The value of R°.(é,7,@) gives us the reflection function of a semi-infinite non-absorbing cloud. It 
has the following simple form at nadir geometry used in this study [8]: 


0.37 + 1.94 cos vo p(n — Vo) 
1 + cos 9o 4(1 + cos 9o)’ 
where p is the cloud phase function. Note, that the spherical albedo r is given by 1 — t for 
non-absorbing clouds and 
— sinh(x + «y)exp(—y) — sinh(y) exp(—x — y) 
sinh(x + ay) 


Roo(Vo) = 


in the presence of absorption. 

The value of g is taken to be equal 0.86, which is typical for clouds. The optical thickness is 
varied in the range 5—100. This is in correspondence with results [3]. The single scattering albedo 
@o is defined by the following equation: 

abs 
Wo = 1—- b 
Oext 


where 
— „A G C 
Fabs = Fabs + Tabs ag Tabs» 
— A G C 
ext = Foxt + O ext a Text? 


where indices A, G and C relate correspondent value to aerosol, gas and cloud extinction or ab- 
sorption coefficients. Let us assume that cf, = oÇ = 0. Then we have 
Mj =1- Fabs 
Cext 
Thus, the value of wp depends on the height. Eq. (A.1) for the cloud reflection is obtained, however, 
for the case of a vertically homogeneous layer. The dependence wo(z) is not particularly strong in 
the area, where cloud is present (see, e.g., the curve wo(z) at the attitude range 4-5 km in Fig. 7). 
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Therefore, we adopt here the model of an “effective homogeneous layer”. In this case one should 
find the height inside the cloud at which the value of wo should be taken for calculations. We 


assume that this height coincides with the Chamberlian’s absorption length [32,33]: 
1 
f= (A.2) 


== $ 
Koext 


where the absorption length / specifies the distance (along the vertical) from the cloud top to a 
fixed level inside the cloud and [32] 


k = y/3(1 — g)( = @o(1/k)). (A.3) 


Solving Eq. (A.3) numerically, we find the value of k, which should be substituted in Eq. (A.2) 
to find the geometrical depth 7, where the single scattering albedo should be taken for calculations 
with Eq. (A.1). We use in Eq. (A.3): g = 0.86. The symbol wo(1/k) means the single scattering 
albedo at the optical depth 1/k from the cloud top. 

Note that Eq. (A.3) is solved by the method of iterations, assuming that @o(1/k) in the first 
iteration is equal to wo at the cloud-top-height. Then Eq. (A.3) is used to find k. This allows to 
find the next approximation for @o(1/k), which is used in Eq. (A.3) to find a new value of k. This 
process is repeated until the convergence is reached. 
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